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Abstract 

This paper reports result of an experimental study in the 
NASA Glenn Icing Research Tunnel (IRT) to evaluate how 
well the current recommended scaling methods developed for 
fixed-wing unprotected surface icing applications might apply 
to representative rotor blades at finite angle of attack. Unlike 
the fixed-wing case, there is no single scaling method that has 
been systematically developed and evaluated for rotorcraft 
icing applications. In the present study, scaling was based on 
the modified Ruff method with scale velocity determined by 
maintaining constant Weber number. Models were unswept 
NACA 0012 wing sections. The reference model had a chord 
of 91.4 cm and scale model had a chord of 35.6 cm. Reference 
tests were conducted with velocities of 76 and 100 kt (39 and 
52 m/s), droplet MVDs of 150 and 195 pm, and with 
stagnation-point freezing fractions of 0.3 and 0.5 at angle of 
attack of 0° and 5°. It was shown that good ice shape scaling 
was achieved for NACA 0012 airfoils with angle of attack up 
to 5°. 

Nomenclature 

A c Accumulation parameter, dimensionless 

b Relative heat factor, dimensionless 

b 0 Stagnation-point relative heat factor 

c Airfoil chord, cm. 

d Cylinder diameter or twice the leading-edge radius of 

airfoil, cm 

h c Convective heat transfer coefficient, cal/sec nr K 

h G Gas-phase mass-transfer coefficient, g/sec m 2 

K Inertia parameter, dimensionless 

K 0 Modified inertia parameter, dimensionless 

L Length proportional to model chord, cm. 

LWC Cloud liquid- water content, g/m 3 

MVD Water droplet median volume diameter, pm 

n Freezing fraction, dimensionless 

w 0 Stagnation-point freezing fraction 

p Pressure, Pa 

p w Vapor pressure of water in atmosphere. Pa 

Pw, Vapor pressure of water at the icing surface. Pa 

Re 5 Reynolds number of water drop, dimensionless 

t Temperature, °C 

tf Freezing temperature, °C 


t s Surface temperature, °C 

T Absolute temperature, K 

V Air velocity, kt 

Weber number based on dimension L and water 

properties, dimensionless 
a Angle of attack, deg. 

P Collection efficiency, dimensionless 

po Stagnation-point collection efficiency 

<(> Droplet energy transfer parameter, °C 

X Drop range, m 

Xstokes Drop range if Stokes Law applies, m 

A f Latent heat of freezing, cal/g 

A v Latent heat of vaporization, cal/g 

ju Air viscosity, poise 

0 Air energy transfer parameter, °C 

p Air density, kg/m 3 

pi Ice density, kg/m 3 

p w Liquid water density, kg/m 3 

a Surface tension of water against air, dyne/cm 

t Accretion time, min 

Subscript 

R Reference 

S Scale 

st Static 

tot Total 

Introduction 

Aircraft and component manufacturers must thoroughly test 
new products to determine the effect of icing on their 
performance. This testing is performed both during the design 
process and for certification purposes. Flight-testing is 
necessary but is expensive and can only be done when 
atmospheric icing conditions exist. Furthermore, it can be very 
time consuming to find in nature the extremes in the cloud 
drop size and liquid-water content envelope required for 
certification testing. Icing wind tunnels can simulate natural 
icing with water-spray and refrigeration systems and provide 
control of cloud conditions, temperature and airspeed to 
permit safe, convenient and relatively inexpensive testing. 
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Because of test-section blockage limitations, many 
components cannot be tested foil size in an icing wind tunnel. 
Furthermore, facilities that simulate natural icing can provide 
only limited ranges of air speed, cloud drop size, and liquid- 
water content. A scaling method is a procedure to determine 
the scaled test conditions to produce the same result as 
exposing the reference model to the desired cloud conditions. 
When the reference (full-size) model is too large for a given 
facility, model-size scaling is applied, and when the desired 
test conditions are outside the facility operating capability, 
test-condition scaling is required. Constraints may also 
prohibit strict geometric scaling of sub-scale models. 

Scaling methods consist of a set of equations that are used 
to determine the necessary scale test conditions, given the 
reference conditions, model size, and geometry that need to be 
simulated. Unlike the fixed-wing case where extensive 
research efforts to develop appropriate scaling methods 
(Refs. 1 and 2) have begun as early as in the 1950s and 
continue to the present, there is no single method that has been 
systematically developed and evaluated for rotorcraft icing 
scalmg applications. The scaling parameters have not been 
adequately verified (Ref. 3), in part because there is no 
publicly available research-quality database with which to 
assess the validity of these methods in a rigorous manner. 
Scaling methodologies also play an important role in obtaining 
data for icing code development and validation (Ref. 4). 

Recently, icing tests were performed in the NASA Glenn 
Icing Research Tunnel (IRT) to begin a preliminary evaluation 
of existing fixed-wing scaling methods for icing applications 
to representative rotor blades at finite angles of attack. All 
tests were made using NACA 0012 models at angles of attack 
of 0° and 5°. The constant 0^ method of determining scale 
velocity was used. Results will be presented for stagnation- 
point freezing fractions of 0.3 and 0.5 and for the two 
reference velocities. In addition, to support an ongoing 
research effort in developing scaling methods for super-cooled 
large droplet (SLD) conditions, it is necessary for the 
reference and scale MVDs to be all in the SLD regime. 


Scaling Equations 

The current recommended scaling methods, from 
References 1 and 2, for fixed-wing unprotected surfaces under 
Appendix C and SLD icing conditions would be a good 
candidate to begin the development of scaling method for 
rotorcraft icing study. These two references provide a detailed 
explanation of the rationale and derivations of the equations 
used to describe the similarity parameters involved in 
traditional fixed-wing aircraft icing. It also includes some 
validation data. A very brief summary of current 
understanding and practice of scaling methods is given here, 
and the reader is referred to these references for a more 
comprehensive discussion. 

For traditional fixed-wing aircraft icing, various icing 
scaling studies over the past 50 years have shown that there 


are four most important similarity parameters affecting ice 
shape. They are: 

(1) Accumulation parameter. 


A = 


LWCVx 
d P,- 


( 1 ) 


(2) Stagnation-point collection efficiency p 0 , for example 
Langmuir and Blodgett (Ref. 5) had shown the p 0 for small 
water drops across a cylinder of radius d to be a function only 
of the modified inertia parameter, K 0 , 


p„=- 


1 - 401^-1 


1+1.401 K„-~ 
0 8 


0.84 


( 2 ) 


This K 0 was defined by Langmuir and Blodgett to describe 
the inertia of drops in an air stream flowing around a body: 


K.-M. 1 


8 X, 


Stokes 


*-! 
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(3) 


In equation (3), K is the drop inertia parameter. 


p MVD 2 V 

_ » u' 

18 d\x 


(4) 


where d is the radius for cylindrical models or twice the 
leading-edge radius for airfoils. Also in equation (3) TJX stokes is 
the drop range parameter, defined as the ratio of actual drop 
range to that if Stokes drag law for solid spheres applied. It is 
a function only of the drop Reynolds number, Re 5 


V MVD p 

Re 5 = 

V 


(5) 


Langmuir and Blodgett tabulated the values of their 
calculated range parameter. For convenience, those data had 
been curve-fitted to the following expression: 


X 

^-Stokes 


0.8388 + 0.001483 Re, 'l 1 

O 

+ 0 . 1847 ^ 


( 6 ) 


(3) Stagnation-point freezing fraction, n 0 . From Messinger’s 
(Ref. 6) surface energy balance, the stagnation-point freezing 
fraction is 


c 

p,ws 


4>+f 

°o J 


(7) 
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The individual terms in this expression are (j), the water energy 
transfer parameter. 


<t> 


_V_ 2_ 
2 c 


0 5 .the air energy transfer parameter. 


0 = 
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P HIV 
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( 8 ) 


( 9 ) 


and 6 0 , the relative heat factor, introduced by Tribus, et 
al.(Ref. 7) At the stagnation line, it is: 


b o = 


LWC VB n c 

r 0 p , 


p.ws 


( 10 ) 


Equation (9) from Ruff (Ref. 8) includes compressibility 
effects. A simpler form without compressibility was used by 
Charpin and Fasso (Ref. 9) and others. Ruffs expression for 0 
was used in the calculations for all this work, but values found 
without compressibility are not significantly different for most 
icing conditions. 

(4) Weber number. 


We, 


™P W 


( 11 ) 


With scale model size selected, by matching scale and 
reference values of We L the scale velocity can be determined. 
By matching p 0 the scale MVD can be found. References 1 and 
2 also showed that the effects of temperature and LWC are not 
independent, but interact through the freezing fraction (i.e., the 
Olsen method). Therefore, with scale LWC chosen, by 
matching n 0 the scale temperature can be calculated. Finally, 
by matching A c the scale accretion time can be established. For 
the scale test, then, only temperature, velocity, MVD and tune 
have to be calculated from the known (reference) values of the 
similarity parameters. 

While some of these similarity parameters are based on 
conditions that apply anywhere on the model, p 0 and n 0 are 
specific to the stagnation lme of a clean model. Therefore, 
strictly speaking, scaling methods only apply at the stagnation 
line of a clean model. These parameters vary with chord-wise 
location and change as ice accretion modifies the geometry. 
Consequently, two assumptions are implied for scaling to be 
valid. The first is that with similar model geometries and 
similar flows around both reference and scale models, if p and 
n match at the stagnation point, they will tend to match 
everywhere on the model. This assumption has been verified 


for collection efficiencies in Reference 2. As for other airflow 
related issues: transition and roughness, for example, may not 
scale, and Re effects are assumed to have a minor influence on 
the final ice shape. Second, if the scaling is done successfully, 
the scale ice shape normalized by the model size will 
consistently agree with the reference for any accretion time 
starting with the clean model. Therefore, scale p and n will 
continue to match the respective reference values, even though 
those parameters are changing with time. 

As for rotorcraft all the similarity considerations for fixed- 
wing aircraft icing should apply. In addition, compressibility 
and dynamic effects due to blade rotation must be considered. 
Also flow-field similarity requires matching of advance ratio 
and pitching settings as well as geometric similarity of the 
rotors. While future studies need to demonstrate scaling when 
aforementioned effects are present, it is necessary for the 
development of the appropriate scaling methods to begin with 
a much simplified flow and model configuration as a baseline 
and then proceed with additional physical effect of importance 
one at a time m testing. This could allow us to evaluate the 
significance of various physical effects unique to rotorcraft 
icing independently and reduces the complications from 
scaling testing. 

Furthermore for rotorcraft icing scaling testing, it may be 
possible to divide the rotor blade into several radial regions 
according to rotating speed, as shown in Figure 1. 
Determination of the local ice shape would come from tests 
performed at various speeds to represent each radial location. 

Therefore in this study the focus will be on the mam rotor to 
address scaling, and the baseline configuration is a 2-D flow 
over a generic rotor blade airfoil at a prescribed angle of attack 
that can simulate local flow conditions along the blade at 
various radial locations. The models chosen are the 
symmetrical NACA 0012 airfoil of different chords. 



Effective 


Main Region 
Of Interest 
(expect rime) 


Figure 1 . — Possible ice accretion regions along the blade. 
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Test Description 

The icing tests were performed in the NASA Glenn Icing 
Research Tumiel (IRT). The IRT is a closed-loop, refrigerated, 
sea-level tunnel with a 1.8 m by 2.7 m rectangular test section. 
The icing cloud is generated by operating 10 spray bars, a 
configuration in use since 1998. 

The IRT cloud calibrations for both Appendix C and SLD 
conditions used for these tests were perfonned in the summer 
of 2008. The TIFC measurements were made using icing blade 
method as reported previously (Ref. 10). However, during the 
2008 calibration there was not enough tune to perform MVD 
measurements. The MVDs reported in this paper are based on 
an analysis of the MVD calibration data completed in 
February, 2006. In addition, because only a few specific 
MVD-LWC combinations at speeds of 100, 150, 200 and 
250 kt (51, 77, 103 and 128 m/s) have been calibrated to date 
in the SLD regime, additional LWC measurements for SLD 
conditions were made in the IRT on September 22, 2008 with 
an icing blade. Therefore SLD tests are constrained to these 
particular conditions. 

The models used were NACA 0012 airfoil sections with 
chords of 91.4 and 35.6 cm. The 91.4-cm-chord airfoil is 
pictured in Figure 2. It was a full-span, fiberglass model and 
served as the reference model. The 35.6-cm-chord scale model 
was of 61 -cm span and made of aluminum. It was mounted 
vertically between splitter plates at the center of the IRT test 
section as shown in Figure 3. The supports permit changing 
the angle of attack. Although no angle-of-attack sweeps are 
planned for any of the present tests, models need to be rotated 
to align with airflow using external pressure belts wrapped 
around the model leading edge at specified chord-wise 
locations along the tunnel vertical center plane. An angle of 
attack of 5° is planned specifically; both reference and scale 
models were rotated to align with the airflow by matching the 
corresponding c p curve using the external pressure taps on the 
models. Horizontal lines at the leading edge indicated tunnel 
center and ±2.5 cm from the center as visual guides for 
locating ice tracings. Both SLD and Mod-1 nozzles were used. 
Also because of the quick start capability of the current IRT 
spray system, the models were not shielded during the 
initiation of the spray. 

In preparing for a test, the temperature and airspeed in the 
test section and the air and water pressures on the spray 
manifolds were set. When these conditions had stabilized, the 
spray nozzle valves were opened to initiate the spray. The 
spray was timed for the required duration, and then turned off. 
The fan was brought to a full stop and the researchers entered 
the test section to document the ice shape with hand tracings. 
Close-up photographs were also taken with a hand-held digital 
camera. 

To record the ice shapes, a thin slice was first melted 
through the ice normal to the model surface. A cardboard 
template was then placed into this slit and an outline of the ice 
shape traced by pencil, giving a two-dimensional cross section 
of the ice. Tracings were taken at the vertical center of the 


tunnel (91 cm from the floor) and at 2.5 cm above the center. 
The ice shapes so recorded were digitized using an automated 
line-following feature in the image-analysis software, 
SigmaScan Pro (Ref. 11). Results from test entries in 2008 and 
2009 will be presented. Since the shape differences between 
the two tracing locations were never significant, only 
centerline shapes will be reported here. 



Figure 2. — 91 .4-cm chord model in the IRT test section. 



Figure 3. — 35.6-cm chord model in the IRT test section. 
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Uncertainty Analysis 

Estimates of the uncertainty in the reported average 
conditions were made by considering inherent errors of 
instruments, temporal fluctuation and spatial variation of the 
instrument readings in the test section, and uncertainty in 
tunnel calibration of MVD and LWC. Recorded air temperature 
was believed to be accurate to ±0.2 °C, although variations 
during the period of an icing spray increases the uncertainty 
for reported average temperatures to about ±0.5 °C. The 
uncertainty in air velocity was estimated to be ±2kt. For 
Appendix-C conditions the net uncertainty in MVD was 
estimated at ±12 percent. For SFD conditions it may have 
been as much as ±20 percent. These uncertainties are not 
referenced to an absolute value of MVD , which is unknown. 
Repeatability and scatter in the LWC calibration data suggests 
the uncertainty is about ±12 percent for both Appendix C and 
SLD conditions. 

The test-parameter uncertainties were used to estimate the 
following uncertainties in the similarity parameters for the 
SLD tests the uncertainties were: 2 percent in p 0 , 12 percent in 
A c , 10 percent in n 0 and 3 percent in 0 L . 

Results 

A preliminary evaluation of the angle of attack effect on 
scaling performance is shown in Figure 4 to Figure 7, which 
compare reference and scale ice shapes obtained for a = 0° 
(part (a) of each figure) with those at a = 5° (part (b) of each 
figure) for stagnation-point freezing fraction of 0.3 and 0.5 at 
reference velocities of 76 and 100 kt. For each figure, 
reference ice shapes are shown shaded, while the scale shapes 
are indicated by a solid line. The table below each figure gives 
the test conditions and similarity parameters for each pair of 
reference and scale tests. The conditions given are the average 
conditions recorded over the duration of each test, which can 
sometimes differ slightly from the planned set points. The 
parameters in the tables were calculated from these average 
conditions. 

Stagnation Freezing Fraction of 0.3 

Figure 4(a) and (b) shows reference and scale ice shape 
comparisons at 0° and 5° AoA with the constant 0 Z method. 
The same reference icing condition was used for both 
comparisons. The reference model size, velocity and MVD 
were 91.4 cm, 76 kt and 147 pm. 

In Figure 4(a) the scale and reference values of p 0 , n 0 , and 0 L 
matched within 5 percent and the values of Po A c were just 
within 8 percent. As demonstrated from previous scaling work 
for a = 0°, the size and shape of the reference ice were well 
simulated by the scale test in the leading-edge region. In Figure 
4(b) the scale and reference values of p 0 , n 0 and 0 L matched 
within 2 percent and the product p oA c agreed only within about 


17 percent. The scale ice shape closely simulated the reference 
main ice shape and feather region, however the size of the scale 
ice was relatively larger than the reference ice. Post data 
analysis indicated that the actual spray time for 09-29-08 run 2 
(see Figure 4) was only 89 percent of the planned value. This 
disagreement explains the difference observed in the leading- 
edge ice thickness for this pair of shapes. 

In addition. Figure 5(a) and (b), gives similar ice shape 
comparison results for a higher reference velocity of 100 kt 
with the 91.4-cm-chord 198-pm-MFD reference conditions 
scaled to 35.6 cm. In both cases (i.e., a = 0° and a = 5°), the 
scale and reference values ofp 0 , P oA c , n 0 and 0 L matched 
within 2 percent. The 35.6-cm-chord scale tests produced 
main ice shapes in fairly good agreement with the reference. 
Even the sizes of large feathers adjacent to the main shape and 
the smaller feathers further aft were simulated reasonably well 
(except for a few larger feathers further aft on the lower 
surface in Figure 5(a)). 

Stagnation Freezing Fraction of 0.5 

Figure 6(a) and (b) each compare reference and scale ice 
shapes with the 91.4-cm reference conditions scaled to 35.6 cm 
at 0° and 5° AoA respectively. The same reference icing 
condition was used for both comparisons. The reference model 
size, velocity and MVD were 91.4 cm, 76 kt and 147 pm. 

In Figure 6(a) the scale and reference values of p 0 , P oA c , n 0 , 
and 0 L all matched within 2 percent. The size and shape of the 
reference ice were well simulated by the scale test in the 
leading-edge region. However, the scale (SLD) feathers were 
larger than the recorded in the SLD reference test. In a 
previous study (Ref. 12) with SLD icing scaling, the larger 
feather formation was observed for the Appendix C scale 
condition. It was speculated then that either there was 
something different about SLD icing physics compared with 
Appendix C or differences in the IRT drop-size distributions 
were responsible. Further study to isolate these large-feather 
conditions is needed to improve scaling performance. In 
Figure 6(b) the scale and reference values of p 0 , n 0 and 0 L 
matched within 2 percent and the product p^ agreed within 6 
percent. The scale ice shape simulated the reference main ice 
shape and feather region fairly good, however the size of the 
scale feathers was still larger and the horn angle of the 
leading-edge main ice shape was also slightly different. 

For a higher reference velocity of 100 kt, similar ice shape 
comparisons were made with the 91.4-cm reference conditions 
scaled to 35.6 cm and the results were shown in Figure 7(a) 
and (b). The reference MVD was 198 pm and the scale MFD’s 
was 87 pm. In both cases, the scale and reference values of p 0 , 
w 0 and We z matched within 2 percent. Reference and scale 
Po^ c in Figure 7(a) matched within about 3 percent, and for 
Figure 7(b) agreed within 1 percent. The scale ice shapes for 
both cases were able to simulate the reference main ice shapes 
and feather regions well. 
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Figure 4. — Scaling from 91 .4 to 35.6-cm-chord with 0 l matched. NACA 0012 airfoils at a = 0 and 5 ; no, 0.3; Vr, 76 kt 
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Figure 5. — Scaling from 91.4 to 35.6-cm-chord with 0 l matched. NACA 0012 airfoils at a = 0° and 5°; no, 0.3; Vr , 100 kt. 
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Figure 7. — Scaling from 91.4 to 35.6-cm-chord with 0 L matched. NACA 0012 airfoils at a = 0° and 5°; n 0 , 0.5; V R , 100 kt. 


NASA/TM — 20 10-215801 


7 

















































































































































These ice shape comparison results would suggest that the 
current recommended scaling methods for fixed-wing icing 
could be used for some generic rotor blade icing applications 
without any modification for some fmite range of angle of 
attack. However, the data is very limited. Further evaluation is 
needed to establish the limiting conditions of such application 
in the IRT to improve scaling performance. 

Conclusion 

Icing tests were performed in the NASA Glenn Icing 
Research Tunnel and scaling results for 0° and 5° AoA were 
presented in this study. Evidence from ice-shape comparison 
shows that for NACA 0012 airfoils the ability to simulate a 
reference ice shape by scaling was not affected by the angle of 
attack in the range tested. Also, for the limited conditions of 
this study, there was no evidence of any difference in the 
fundamental formation mechanisms of either feather growth or 
horn formation on airfoil models at angle of attack. 

Good ice shape scaling was achieved in this study by 
matching scale and reference values of the parameters p 0 , A c , 
«o and 0 L for NACA 0012 airfoils with angle of attack up to 
5°. Model size ratio was 2.6:1 and freezing fractions covered 
the range from 0.3 to 0.5. The present SLD reference tests 
were made with velocities of 76 and 100 kt, and these 
conclusions may not be valid for higher velocities. Similar 
concern for larger static angle of attack should also be 
considered. Additional testing in the IRT is needed to 
document the limiting conditions and to evaluate other 
possible alternatives of scaling methods for such application to 
improve the current scaling knowledge and capability for 
rotorcraft icing. 

Further icing testing on oscillating airfoils is planned in the 
IRT in coming years. This testing should provide a better 
evaluation of scaling performance in a much more realistic 
rotor blade operation environment. 
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